Quantum Reservoir Engineering 



J.F. Poyatos*, J.I. Cirac,* and P. Zoller 

Institut fur Theoretische Physik, Universitat Innsbruck, Technikerstrasse 25, A-6020 Innsbruck, 

Austria. 
(February 5, 2008) 

Abstract 

We show how to design different couplings between a single ion trapped in a 
harmonic potential and an environment. This will provide the basis for the 
experimental study of the process of decoherence in a quantum system. The 
coupling is due to the absorption of a laser photon and subsequent sponta- 
neous emission. The variation of the laser frequencies and intensities allows 
one to "engineer" the coupling and select the master equation describing the 
motion of the ion. 



One of the fundamental questions of Physics is understanding the borderline between 
microscopic phenomena ruled by quantum mechanics and the macroscopic world of classical 
physics. In particular, according to quantum mechanics [|1]], a system can exist in a superpo- 
sition of distinct states, whereas these superpositions seem not to appear in the macroscopic 
world. One possible explanation of this paradox is based on the fact that systems are 
never completely isolated but interact with the surrounding environment, that contains a 
large number of degrees of freedom. The environment influences the system evolution which 
continuously decoheres and transforms system superpositions into statistical mixtures which 
behave classically fHH. This problem is directly related to the problem of measurement in 
quantum theory 0|| where the system to be measured is described by quantum mechanics 
and the measurement apparatus is assumed to behave classically. Apart from this funda- 
mental point of view a more practical aspect is the question to what extent one can preserve 
quantum superpositions, which is the basis of potential applications of quantum mechanics, 
such as quantum cryptography and computation 

From the theoretical point of view, quantum decoherence has been studied extensively 
PJ3|J^[TT|. Most of the effort has been focused on the decoherence of a harmonic oscillator 
(the system) due to the coupling to a reservoir consisting of oscillators (the environment). 
According to these studies, the decoherence process depends critically on the form of the 
coupling between the system and the environment. On the experimental side, however, there 
have not been any systematic investigations. This is due to the lack of experimentally acces- 
sible systems where both the time scale and the form of the system-environment coupling 
can be changed in a controlled way. 

In this paper we will show how to "engineer" the system-environment coupling in a 
situation that is experimentally accessible with existing technology. The system of interest 
will be an ion confined in a electromagnetic trap, and the environment will be the vacuum 
modes of the electromagnetic field. This corresponds to an experimental realization of a 
harmonic oscillator coupled to a reservoir of oscillators. The coupling between our system 
and the environment takes place through the recoil experienced by the ion when it inter- 
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changes photons with the electromagnetic field. As we will show below, this coupling can be 
manipulated by laser radiation. Variations of the laser frequency and intensity allows one 
to engineer such a coupling. 



Laser cooled trapped ions |12| are a unique experimental system: unwanted dissipation 
can be made negligible for very long times, much longer than typical times in which an 
experiment takes place. Furthermore, arbitrary quantum states of the ion's motion can be 
synthesized and coherently manipulated using laser radiation fl3[| . In addition, the state of 



motion can be completely determined in the sense of tomographic measurements ||14j| . In a 



series of remarkable experiments, Wineland and collaborators have generated a variety of 



non-classical states of ion motion [15,19]. In particular, they have been able to produce [16 



a so-called "Schrodinger cat state" [17J corresponding to 



oc | a) + | - a), (1) 

with \a) = J2T=o ce n /Vn\\n) a coherent (quasiclassical) state. In Fig. 1(a) we have plotted 
the density operator for such a state in the position representation, i.e. (the real part of) 
p(x,x') = (x\p\x'). The peaks near the diagonal correspond to two possible localizations of 
the particle, whereas the other two peaks are related to the coherences that are responsible 
from the quantum behavior 0. 

The interaction of a Schrodinger cat state with the environment has been the paradigm 
of decoherence of superposition states. As first argued by Zurek |2| (see also Refs. |||2|||U), 
for a coupling which is linear in the system coordinates, a macroscopic superposition of the 
form (Q) decays to a statistical mixture p oc + | — a)(—a\, on a short time scale 

(decoherence time) which is related to the size of the cat (|a| 2 ) and is much faster than the 
energy dissipation time: this provides an explanation for the absence of superpositions in the 
macroscopic world ||. We emphasize that the decoherence process of ([I]) is sensitive to the 
form of the reservoir coupling. For some quadratic couplings, for example, the decoherence 
and energy dissipation time can become identical |S]]; moreover, there exist interactions 
which allow Schrodinger cat states to be stable, and, what is more surprising, dissipation 



can drive a system into a steady state of the form ([!]) @. For example, in Fig. l(b,c) the 
decay of a Schrodinger cat under linear and quadratic coupling is illustrated: for a linear 
coupling (Fig. 1 b) the nondiagonal peaks (coherences) of the density matrix decay much 
faster than for the quadratic couplings (Fig. 1 c). We will show that all these theoretical 
predictions can be tested experimentally for the case of a trapped ion. 

The process of decoherence can be analyzed in detail under very general assumptions 
invoking to the so-called Markov approximation, which considers the correlation time for the 
environment to be much shorter that the evolution time of the system due to the coupling 
| TT| . In this case the interaction of a system with an environment is described in terms of a 
master equation. For a single decoherence channel this equation has form (H = 1) 

p = 7 (2/p/t-/t / p- p/ t /)) (2 ) 

Here p is the reduced density operator for the system in the interaction picture after tracing 
over the reservoir. The operator / and the parameter 7 reflects the system-environment 
coupling. For a harmonic oscillator / will be a function of the creation and annihilation 
operators a and a*, which are defined as usual X = l/{2Mu) 1 ' 2 (a^+a), P = i(Mv/2) 1 ' 2 (a ! — 
a) where X and P are the position and momentum operators and M the particle's mass. 
According to Zurek , the coupling with the environment singles out in a quantum system 
a preferred set of states, sometimes called "the pointer basis" . This basis depends on the 
form of the coupling /. For example, for / = X the pointer basis is the position eigenstates. 
The density operator describing the system evolves in such a way that it rapidly becomes 
diagonal in this preferred basis, which is usually connected to the disappearance of quantum 
interferences. Our goal is now to find an experimental realization of the master equation (0) 
for different system-reservoir couplings / = f(a,a'). 

Let us consider a single ion moving in a one-dimensional harmonic potential. The ion 
interacts with a laser in a standing wave configuration of frequency 00 l, close to the tran- 
sition frequency u of two internal levels \g) and |e). Using standard methods in quantum 
optics based on the the dipole, Born-Markov, and rotating wave approximations, the master 
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equation that describes this situation can be written in the general form 

p = -iH cS p + ipHl s + Jp, (3) 

with 

r 

H eff = H tp + H int + H cou - i-|e)(e|, (4a) 
Jp = T J 1 duN{u)e- iriu{a+a1) (j-p(j + e iriu{a+a]) (4b) 

where H tp = ua^a, H int = \u a z , and H cou = f sin[?7(a + a f ) + (f)](a + e' luJLt + <j_e luJLt ), give 
the free Hamiltonian for the motion in the trap, the internal two-level system Hamiltonian, 
and the one describing the coupling with the lasers, respectively. Here, a + = \e)(g\ = (&-)' 
and a z = \e)(e\ — \g){g\ are usual spin-^ operators describing the internal transition, v is the 
trap frequency, T the spontaneous emission rate, and r] = (/c| / /2Mz/)( 1 / 2 ^ the Lamb-Dicke 
parameter which is the ratio of the recoil frequency k 2 L /2M to the trap oscillation frequency 
v. In the expression for the superoperator J , the exponentials are related to the photon 
recoil that takes place in each spontaneous emission process, and the integral takes into 
account the different angles at which that photon can be emitted, with a normalized dipole 
pattern N(u). In the Hamiltonian describing the coupling with the lasers, Q is the Rabi 
frequency, and <fi characterizes the relative position of the trap center with respect to the 
node of the laser standing wave. Here we will assume that either = (excitation at the 
node of the standing wave) or = tt/2 (excitation at the antinode). 

We will proceed now by simplifying the master equation for the ion in a regime defined by 
three limits which are typically fulfilled in experiments |I5| , |16[ : (i) Lamb-Dicke, (ii) strong 
confinement (iii) low intensity. The first one allows to expand the above master equation in 
terms of the Lamb-Dicke parameter 77 <C 1, retaining only the orders that contribute to the 
dynamics. The second one assumes r v and together with the third one allows to include 
in the coupling Hamiltonian only on-resonance terms (secular approximation). Finally, the 
third one assumes a sufficiently low laser intensity (the specific form of this limit will be 
given later), and will serve us to adiabatically eliminate the internal excited level |e). 
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Let us start by simplifying the coupling Hamiltonian under the above limits. To do that, 
we move to a rotating frame defined by the unitary operator U = e~ i ^ Htp+Hint ' t . Following 



Ref. |13[ we assume that: (i) For excitation at the node (0 = 0), 5 = ul — = (2/c + 1)^ 
(k = 0, ±1, . . .) (ii) For excitation at the antinode {(f) = vr/2), 5 = 2kv (k = 0, ±1, . . .). In 
this rotating frame, after performing the rotating wave approximation and the Lamb-Dicke 
expansion, we obtain H cou = ^ (<?"+/ + Pc-)) where both Q' and the form of the operator 
/ depend on the frequency of the laser. For example, for 5 = —u, we have f = a, and 
Q' = Qrj/2, whereas for S = —2u, f = a 2 , and Q' = —Qi] 2 /6. Apart from the strong 
confinement, in the first case, the secular approximation can be performed for Q' u, 
whereas in the second case it is needed Q 2 /u Q'. This two conditions can always be 
fulfilled for low enough laser intensity, and together with fi' < T define the low intensity 
limit. 

In the next step we eliminate adiabatically the internal excited state using standard 
procedures of quantum optics [LT]|. Physically, since H' < T the ion practically spends no 



time in the excited level and therefore we can eliminate it. Finally, expanding in powers 
of 7] we find the desired master equation (|^), with corrections of the order rf. The master 
equation will be valid for times such that these corrections are not important, that is for times 
t ( / yr] 2 n)~ 1 where n is the typical phonon number of the state of the ion. Nevertheless 
in the Lamb-Dicke limit this time can be much longer than the time required to reach 
the steady state using the approximated master equation. Note that through the adiabatic 
elimination we are coupling effectively the motion of the ion with the environment. The 
fact that this coupling takes place through the absorption of laser photons, and we have 
to choose the way on how this actually happens, allows one to manipulate the coupling 
system-environment . 

According to our analysis, by varying the laser frequency we obtain the master equation 
(0) with different coupling operators /. In Fig. 2 we have illustrated the laser configurations 
which produce several / operators. In Fig. 2(a), for example, the laser is tuned to the so- 
called "lower motional sideband," 5 = —u, and the ion is located at the node of the standing 



wave field which leads to a coupling operator f = a. This can be easily understood by 
noting that in each absorption and spontaneous emission cycle one phonon is annihilated 
on a time scale given by the optical pumping time. Similarly, in Fig. 2(b) the laser is 
tuned to the "second lower sideband" 5 = —2u at the antinode of the laser standing wave 
which gives the two— phonon coupling f = a 2 . These two cases of linear and quadratic 
coupling correspond to the two examples discussed in Figs. l(b,c). In fact, these figures 
were obtained by a numerical solution of the full master equation @ with quantum Monte- 
Carlo wavefunctions simulations [|18|]. As noted before, the decoherence acts in a different 
way depending on the coupling operator, according to our previous discussion. 

It is simple to generalize the above derivation to find situations with other interesting 
(and perhaps unusual) coupling operators /. For example, consider the case in which two 
lasers of frequency u + v and Uq — v interact with the ion [Fig 2(c)]. This corresponds to a 
coherent excitation of the lower and upper motional sidebands [|3|]. In this case, following the 
same arguments, one can easily show that the operator is / = fia+ua^ , where fi 2 — u 2 = 1 and 
n/v is the quotient of the Rabi frequencies. This operator corresponds to a squeezed vacuum 
coupling which has been the basis for numerous theoretical predictions in quantum optics 
]TT[ . In particular, choosing equal Rabi frequencies, the coupling is / = a + oc X. This 
corresponds to the case analyzed theoretically by Unruh and Zurek, Caldeira and Legget, 
and other authors p],|3| to describe the decoherence process in terms of the projection of the 
state of the system onto the pointer basis given, in this case, by the position eigenstates. 
Another interesting combination of lasers [Fig. 2(d)] yields / = (a — a)(a — j3), where a and 
(3 are given complex numbers. For a = —f3 the Schrodinger cat state ([I]) is an eigenstate 
of this operator with zero eigenvalue, and thus this state does not decohere under this form 
of coupling. As an aside we note that one can employ this particular form of system- 
reservoir coupling to generate a cat state ([I]) by choosing as the initial state the ground level 
|0) PJ11|. Tuning a laser on resonance at the antinode of a standing light wave one can 
design the coupling in the form of a quantum nondemonolition measurement of the phonon 
number, / = a^a, with the Fock states as the pointer basis. Using more complicated laser 



configurations one can readily show that, for example, /-operators like a(a)a — n) can be 
engineered. The operator a) a — n projects out a subspace with exactly n phonons. Finally, 
exotic combination of operators such as a(a^a — n) — ^(a^a — n — 1)) — 1 can be realized, 
such that the system is driven into the superpositon state \n) + \n+ 1). 

Obviously, there are numerous possibilities to generalize the concept of reservoir engi- 
neering in ion traps. First of all, decoherence of a two or three-mode system can be studied 
by considering the two or three dimensional motion of a trapped ion, respectively. Fur- 
thermore, a master equation with more than one decoherence channel, i.e., an equation 
containing sums of damping terms of the form (0) with different operators fa (i — 1, . . . , N) 



11], can also be easily implemented . This can be accomplished by exciting transitions 



with several incoherent lasers. Another important generalization concerns the possibility of 
coupling a two-level system to a harmonic oscillator (Jaynes-Cummings model) which in 
turn is coupled to an environement. In particular, this will allow to test experimentally one 
of the outstanding predictions of quantum optics [p0[] , namely the damping of a two-level 
system interacting with a squeezed reservoir. Finally, these ideas can be extended to linear 
ion traps in order to study collective effects in an iV-atom + harmonic oscillator system. 

In summary, we have shown how the coupling of a harmonic oscillator (represented by 
the motion of a trapped ion) to an environment can be engineered. We believe that this 
opens a new field in the sense that it will allow for the first time to study experimentally in 
a controlled and systematic way the effects of decoherence in a quantum system. 
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FIGURES 

FIG. 1. (a) p(x, x 1 ) for the state (|]) (a = 3). (b,c) Numerical simulation of the interaction with 
a laser for a time r = 0.067~ 1 and rj = 0.03: (b) ujl = ujq — u, (/ = a); (c) u>l = ujq — 2u, (/ = a 2 ). 

FIG. 2. Laser configurations for several coupling operators /. (a) Laser tuned to 
|n, g) — > \n — l,e), which rapidly decays into the state \n — l,g) leading to / = a. (b) / = a 2 , (c) 
f = fia + va) . (d) / = (a — a) (a — 
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This figure "figl.gif" is available in "gif" format from: 



http://arXiv.org/ps/atom-ph/9603002vl 



